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Abstract 

The relativistic mean-field theory provides a framework in which 
the nuclear many-body problem is described as a self-consistent sys- 
tem of nucleons and mesons. In the mean-field approximation, the 
self-consistent time evolution of the nuclear system describes the dy- 
namics of collective motion: nuclear compressibility from monopole 
resonances, regular and chaotic dynamics of isoscalar and isovector 
collective vibrations. 



1 The relativistic mean field model 

Relativistic mean-field (RMF) models have been successfully applied in cal- 
culations of nuclear matter and properties of finite nuclei throughout the pe- 
riodic table. In the self-consistent mean-field approximation, detailed calcu- 
lations have been performed for a variety of nuclear structure phenomena [|I|]. 
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In this work we present applications of RMF to the dynamics of collective vi- 
brations in spherical nuclei. In relativistic quantum hadrodynamics 0, the 
nucleus is described as a system of Dirac nucleons which interact through 
the exchange of virtual mesons and photons. The Lagrangian density of the 
model is 

C = ${ii-d-m)il> + ^{da) 2 -U{a) 
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The Dirac spinor ip denotes the nucleon with mass m. m a , m u , and m p are 
the masses of the cr-meson, the tu- meson, and the p- meson, and g a , g u , and 
g p are the corresponding coupling constants for the mesons to the nucleon. 
U(a) denotes the nonlinear a self-interaction, and Q^ u , R^ u , and F^ u are 
field tensors [Q, Qj. The coupled equations of motion are derived from the 
Lagrangian density (H). The Dirac equation for the nucleons: 

id t ipi = 
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and the Klein-Gordon equations for the mesons: 

(d 2 - A + ml) a = -g a p s - g 2 o 2 - g 3 a 3 (3) 

(dt-A + mfyur = g u]p (4) 

(dt-^ + ™ 2 p )p, = g p % (5) 

(d 2 t -A)A p = ef™. (6) 

In the relativistic mean-field approximation, the nucleons described by single- 
particle spinors ipi are assumed to form the A-particle Slater determinant |$), 
and to move independently in the classical meson fields. The sources of the 
fields, i.e. densities and currents, are calculated in the no-sea approxima- 
tion B: the scalar density: p s = Y.t=i i'i'^ii ^he isoscalar baryon current: 
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j M = Yld^^il^^i-, the isovector baryon current: = Y^x^il^^i-, 
the electromagnetic current for the photon-field: j£ m = T^x^il^^-f^^i- 
The summation is over all occupied states in the Slater determinant |$). 
Negative-energy states do not contribute to the densities in the no-sea ap- 
proximation of the stationary solutions. It is assumed that nucleon single- 
particle states do not mix isospin. 

The ground state of a nucleus is described by the stationary self-consistent 
solution of the coupled system of equations for a given number of 

nucleons and a set of coupling constants and masses. The solution for the 
ground state specifies part of the initial conditions for the time-dependent 
problem. The dynamics of nuclear collective motion is analyzed in the frame- 
work of time-dependent relativistic mean-field model, which represents a rel- 
ativistic generalization of the time-dependent Hartree-Fock approach. For a 
given set of initial conditions, i.e. initial values for the densities and currents, 
nuclear dynamics is described by the simultaneous evolution of A single- par- 
ticle Dirac spinors in the time-dependent mean fields. The equations 
are equivalent to the equation of motion for the one-body density operator 

P = P(t) 

d 

ift-Q7P=[hD,p\, (7) 

with an initial condition for p 

f>{Un) = Pin- (8) 

ho is the single-nucleon Dirac hamiltonian defined in Eq. @. Starting from 
the self-consistent solution that describes the ground-state of the nuclear 
system, initial conditions are defined to simulate excitations of giant reso- 
nances in experiments with electromagnetic or hadron probes. Frequencies of 
eigenmodes are determined from a Fourier analysis of dynamical quantities. 
In this microscopic model, self-consistent time-dependent mean-field calcu- 
lations are performed for multipole excitations. An advantage of the time- 
dependent approach is that no assumption about the nature of a particular 
mode of vibrations has to be made. Retardation effects for the meson fields 
are not included in the model, i.e. the time derivatives <9 t 2 in the equations 
of motions for the meson fields are neglected. This is justified by the large 
masses in the meson propagators causing a short range of the corresponding 
meson exchange forces. Negative energy contributions are included implicitly 
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in the time-dependent calculation, since the Dirac equation is solved at each 
step in time for a different basis set 0. Negative energy components with re- 
spect to the original ground-state basis are taken into account automatically, 
even if at each time step the no-sea approximation is applied. 

The description of nuclear dynamics as a time- dependent initial-value 
problem is intrinsically semi-classical, since there is no systematic procedure 
to derive the initial conditions that characterize the motion of a specific mode 
of the nuclear system. The theory can be quantized by the requirement that 
there exist time-periodic solutions of the equations of motion, which give 
integer multiples of Planck's constant for the classical action along one pe- 
riod [|J. For giant resonances the time-dependence of collective dynamical 
quantities is actually not periodic, since generally giant resonances are not 
stationary states of the mean-field Hamiltonian. The coupling of the mean- 
field to the particle continuum allows for the decay of giant resonances by 
direct escape of particles. In the limit of small amplitude oscillations, how- 
ever, the energy obtained from the frequency of the oscillation coincides with 
the excitation energy of the collective state. In Refs. |3], [|, ^] we have shown 
that the model reproduces experimental data on giant resonances in spherical 
nuclei. 

2 Monopole resonances and nuclear compress- 
ibility 

The study of isoscalar monopole resonances in nuclei provides an important 
source of information on the nuclear matter compression modulus K nm . This 
quantity is crucial in the description of properties of nuclei, supernovae ex- 
plosions, neutron stars, and heavy-ion collisions. In principle the value of 
K nm can be extracted from experimental energies of isoscalar monopole vi- 
brations in nuclei (giant monopole resonances GMR). However, the complete 
experimental data set on isoscalar GMR does not limit the range of K nm 
to better than 200 — 350 MeV. Microscopic calculations of GMR excitation 
energies might provide a more reliable approach to the determination of the 
nuclear matter compression modulus. Modern non- relativist ic Hartree-Fock 
plus RPA calculations, using both Skyrme and Gogny effective interactions, 
indicate that the value of K nm should be in the range 210-220 MeV. In rela- 
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tivistic calculations on the other hand, both time-dependent and constrained 
RMF results indicate that empirical GMR energies are best reproduced by 
an effective force with K nm ^ 250 — 270 MeV. 

In Ref. || we have performed time- dependent and constrained RMF cal- 
culations for monopole giant resonances for a number of spherical closed shell 
nuclei, from 16 to the heavy nucleus 208 Pb. For the effective Lagrangian we 
have used six parameter sets, which differ mostly by their prediction for 
K nm , but otherwise reproduce reasonably well experimental data on nuclear 
properties. The idea is to restrict the possible values of the nuclear mat- 
ter compression modulus, on the basis of the excitation energies of giant 
monopole states calculated with different effective interactions. In addition 
to the four non-linear sets NL1, NL3, NL-SH and NL2, we have also in- 
cluded two older linear parametrizations, HS and LI. The sets NL1, NL-SH 
and NL2 have been extensively used in the description of properties of finite 
nuclei In order to bridge the gap between NL1 (K nm = 211.7 MeV), 
and NL-SH (K nm = 355.0 MeV), we have also included a new effective in- 
teraction NL3 (K nm = 271.8 MeV). This new parameter set provides 
an excellent description not only for the properties of stable nuclei, but also 
for those far from the valley of beta stability. From the energy spectra and 
transition densities calculated with these effective forces, it has been possi- 
ble to study the connection between the incompressibility of nuclear matter 
and the breathing mode energy of spherical nuclei. For the isoscalar mode 
we have found an almost linear relation between the excitation energy of 
the monopole resonance and the nuclear matter compression modulus. For 
the determination of K nm especially relevant are microscopic calculations of 
GMR excitation energies in heavy nuclei. The results of TD RMF calcula- 
tions for 208 Pb are displayed in Fig. 1: time- dependent monopole moments 
(r 2 (t)) = ^ ($ (t) | r 2 1$(£)) and the corresponding Fourier power spectra for 
the nonlinear effective interactions. As one would expect for a heavy nucleus, 
there is very little spectral fragmentation and a single mode dominates, at 
least for NL1 and NL3. The experimental excitation energy 13.7 ± 0.3 MeV 
is very close to the frequency of oscillations obtained with the NL3 parame- 
ter set: 14.1 MeV. The calculated excitation energy for the NL1 parameter 
set (-ftr nm = 211.7 MeV), is approximately 1 MeV lower than the average 
experimental value. For the linear effective forces HS and LI the oscillations 
are more anharmonic, and the monopole strength is located well above the 
experimental GMR energy. 
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The effective interactions NL1 and NL3 seem to produce GMR excitation 
energies which are quite close to the experimental values. For these two pa- 
rameter sets we have calculated the isoscalar giant monopole resonances in 
a number of closed-shell nuclei: 40 Ca, 56 Ni, 100,114,132^ 90,122^ i46 Gd The 
results are shown in Fig. 2. The energies of giant monopole states are deter- 
mined from the Fourier spectra of the time-dependent monopole moments, 
and are displayed as function of the mass number. The NL1 excitation en- 
ergies are systematically lower, but otherwise the two effective interactions 
produce very similar dependence on the mass number. The empirical curve 
E x « 80 A~ 1/3 MeV is also included in the fi gure, and it follows very closely 
the excitation energies calculated with the NL3 parameter set. Similar results 
are obtained from constrained RMF calculations. Both methods indicate 
that, in the framework of relativistic mean field theory, the nuclear mat- 
ter compression modulus K mn 250 — 270 MeV is in reasonable agreement 
with the available data on spherical nuclei. This value is approximately 20% 
larger than the values deduced from recent non-relativistic density dependent 
Hartree-Fock calculations with Skyrme or Gogny forces. It should be also 
emphasized that the excitation energy of the isovector monopole resonance 
in 208 Pb, calculated with the NL3 effective force, is in excellent agreement 
with the experimental value 26 ± 3 MeV. 



3 Nonlinear dynamics of giant resonances 

Atomic nuclei provide excellent examples of quantum systems in which the 
transition from regular to chaotic dynamics can be studied. Signatures of 
chaotic dynamics have been observed in correlations of nuclear level distribu- 
tions, and in the microscopic and collective motion of the nuclear many-body 
system. Especially interesting in this respect are giant resonances: highly col- 
lective nuclear excitations whose properties, excitation energies and widths, 
nevertheless reflect the underlying microscopic dynamics. Theoretical studies 
have shown that regular collective modes coexists with chaotic single-nucleon 
motion: the adiabatic mean-field created by the nucleons averages out the 
random components of their motion. However, it has also been shown that 
important differences exists between the isoscalar and isovector collective 
modes. 

In particular, we have studied isoscalar and isovector monopole oscilla- 
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tions in spherical nuclei M , but analogous considerations apply to higher mul- 
tipolarities. In Fig. 3 results are shown of time- dependent relativistic mean- 
field calculations for isoscalar and isovector oscillations in 208 Pb. Calculations 
have been performed for the NL1 effective interaction. In the isoscalar case 
both proton and neutron densities are radially expanded, while for the isovec- 
tor mode the proton density is initially compressed by the same amount. In 
Fig. 3a we plot the time series of the isoscalar monopole moment (r 2 (t)), 
and in Fig. 3b the corresponding isovector moment < r^(t) > — < r 2 (t) > 
is displayed. The isoscalar mode displays regular undamped oscillations, 
while for the isovector mode strongly damped anharmonic oscillations are 
observed. On the right-hand panels we plot the corresponding Fourier power 
spectra. The Fourier spectrum of the isovector mode is strongly fragmented, 
but the main peaks are found in the energy region 25 — 30 MeV, in agreement 
with the experimental data. For the isoscalar mode, the time series of the 
monopole moment and the Fourier spectrum show that the oscillations of 
the collective coordinate are regular. On the other hand, the appearance of 
a broad spectrum of frequencies seems to indicate that the isovector oscilla- 
tions are chaotic. For time-series that result from linear physical processes 
the Fourier analysis unfolds the characteristic frequencies which are invari- 
ants of the dynamics, i.e. they classify the dynamics. For nonlinear systems 
the corresponding analysis is somewhat more complicated. 

We have performed the analysis of phase spaces reconstructed from time- 
series of collective dynamical variables that characterize the isoscalar and 
isovector oscillations. For the reconstruction of the dynamics two principal 
quantities have to be determined: the time delay and the dimension of the 
phase space on which the attractor unfolds. The delays which define the 
time-lagged variables have been determined from the first minima of the av- 
erage mutual information function: 27 fm/c for the isoscalar, and 13 fm/c 
for the isovector mode. For the embedding dimensions the method of false 
nearest neighbors has been used: d^ — 3 (isoscalar mode), and d^ — 4 
(isovector mode). The reconstructed phase space can be represented by a 
recurrence plot. For the isoscalar mode the recurrence plot displays a pat- 
tern representative for regular oscillations, for the isovector mode it indicates 
non-stationarity ||. If the dynamics of a system is deterministic, the ensem- 
ble of phase space trajectories converges towards an invariant subset of the 
phase space - the attractor. For chaotic dynamics the attractor has fractional 
dimension, whereas the dimension is integer for regular dynamics. The cor- 
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relation dimension D 2 of the attractor can be numerically evaluated from the 
log- log plot of the correlation integral Ci{r) vs the distance r in phase space, 
for a set of increasing dimensions of the phase space. As the embedding 
dimension increases, D 2 should saturate at a value equal to the attractor's 
correlation dimension (Fig. 4). For the isoscalar mode, for d > 3, the corre- 
lation dimension saturates at D 2 = 2; it does not saturate for the isovector 
mode, but slowly increases to some fractional value between 2 and 3. The 
fractional dimension of the attractor would imply chaotic or stochastic dy- 
namics. 

The identification and quantification of the regular or chaotic dynamics 
can be also based on the evaluation of information-theoretic functionals. For 
the time-dependent one-body nucleon densities, we have calculated the von 
Neumann information entropy functionals. The Fourier analysis has shown 
that the entropy of the isoscalar mode contains the same information as the 
dynamical variable, but the structure is more complicated for the isovector 
mode. Giant resonances represent collective oscillations of densities with 
finite spatial extension, and therefore provide excellent physical examples 
for the analysis of systems that have spatial as well as temporal structure. 
For a nonlinear system in chaotic regime, the influence of spatial motion on 
temporal chaos cab be studied: what are the spatial correlations in a finite 
system that displays chaotic oscillations of a collective dynamical variable? 
These correlations have been described with a time-dependent conditional 
entropy defined from a two-body nucleon density, and the result is a high 
degree of two-body correlations for the isoscalar mode, and very little spa- 
tial correlation for the isovector oscillations of the nucleon density || . From 
the dynamical variables that characterize the proton and neutron distribu- 
tions, we have calculated the average mutual information for the isoscalar 
and isovector modes (Fig. 5). This function quantifies the information that 
is contained in the dynamical variable of the neutron distribution, about 
the proton subsystem, and vice versa. The two curves that correspond to 
isoscalar and isovector oscillations, are plotted as functions of the size of the 
box in the linear embedding of the time-series. The average information that 
is contained in the collective dynamical variable of the proton density, about 
the neutron density, is more than a factor three larger for the isoscalar mode. 

Time-series, Fourier spectra, phase-space plots, Poincare sections, auto- 
correlation functions, and Lyapunov exponents |7| have also shown that the 
motion of the collective coordinate is regular for isoscalar oscillations, and 
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that it becomes chaotic when initial conditions correspond to the isovector 
mode. 
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Figure Captions 



Fig. 1 Time-dependent isoscalar monopole moments < r 2 > (t) and the 
corresponding Fourier power spectra for 208 Pb. The parameter sets are 
NL1, NL3, NL-SH and NL2. 

Fig. 2 Excitation energies of isoscalar giant monopole resonances in spher- 
ical nuclei as function of the mass number. The effective interactions 
are: NL1 (squares) and NL3 (circles). The solid curve corresponds to 
the empirical relation ps 80 A' 1 ^ MeV. 

Fig. 3 Time-dependent isoscalar and isovector monopole moments for 208 Pb. 

Fig. 4 Correlation dimension D 2 as function of embedding dimension, for 
isoscalar and isovector oscillations. 

Fig. 5 Mutual information between the time- dependent mean square radii 
of the proton and neutron density distributions. 
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